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Abstract 

We investigate a disordered version of a thermodynamic fiber bundle model proposed by Selinger, 
Wang, Gelbart, and Ben-Shaul a few years ago. For simple forms of disorder, the model is ana- 
lytically tractable and displays some new features. At either constant stress or constant strain, 
there is a non monotonic increase of the fraction of broken fibers as a function of temperature. 
Moreover, the same values of some macroscopic quantities as stress and strain may correspond to 
different microscopic configurations, which can be essential for determining the thermal activation 
time of the fracture. We argue that different microscopic states may be characterized by an exper- 
imentally accessible analog of the Edwards- Anderson parameter. At zero temperature, we recover 
the behavior of the irreversible fiber bundle model. 

PACS numbers: 62.20.Mk, 64.60.My, 05.50.+q 



I. INTRODUCTION 



It is known that some features of the phenomena of fracture can still be captured by 
simple and schematic models . The simple Fiber Bundle (FB) model, introduced by 

Peirce j4| and Daniels to account for the breakdown load in textiles, is a simple model 
that has provided a useful general scheme for describing the fracture of disordered systems 
in the quasi-static regime. It indicates the critical onset [6] of the fracturing process 0,0], 
and shares a number of features with other models, including a correspondence with results 
from a mean-field approach [9| that is very effective to deal with elasticity 

Macroscopic fracture is a global irreversible process. At the microscopic level, however, 
microcracks can heal under suitable conditions, which gives room to the consideration of 
reversible or at least partially reversible models. In delayed (activated) fractures, the failure 
load may decrease as we increase the temperature or the time of loading . Delayed frac- 
ture phenomena have been investigated by using both fully and partially reversible models 



11 



12 



13 



14 



15 



0, IZ.Q, Q. 



In this article, we study the behavior of an extension 



of a reversible thermodynamical model that has been proposed a few years ago by Selinger, 
Wang, Gelbart, and Ben-Shaul {2^. The model consists of a bundle of independent and 
identical fibers (elastically harmonic strings, with the addition of a failure energy level), 
which are loaded in parallel, and whose states, broken or intact, are allowed to fluctuate 
according to (reversible) Boltzmann weights. Despite its enormous simplicity, this model 
still captures some relevant and essential features of homogeneous nucleation of fractures in 
defect-free crystals, with relevant predictions for systems as iron whiskers 2l(| , which remain 
intact, in a metastable state, until eventually undergoing a fracture. 

It has been recently pointed out that delayed fracture can also occur in non homogeneous 
materials [l|, [lol, [3]. We then decided to investigate the Thermodynamic Fiber Bundle 
(TFB) model of Selinger and coworkers in the presence of disorder. Although we recognize 
that the full description of failure requires the introduction of some degree of irreversibility, 
even at finite temperatures, we believe that the present investigation is a first step towards an 
understanding of the underlying processes of thermally activated fractures. In several sim- 
ple cases, the general thermodynamic properties of this Disordered Thermodynamic Fiber 
Bundle (DTFB) model are analytically accessible, and may differ from the corresponding 
results for the FB and the homogeneous TFB models. In particular, depending on the type 
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of disorder, the fraction of intact fibers can display a non-monotonic behavior with temper- 
ature and constant apphed strain (or stress). The same macroscopic load (strain, stress) 
thus corresponds to different microscopic configurations, and this may have an important 
influence on the activation time of the fracture. We conjecture that these configurations 
may be distinguished by an analog of the Edwards- Anderson parameter, and that the value 
of this parameter can be measured from macroscopic quantities. 

In Sec. II we introduce the TFB model, and recover the original results of Selinger and 
coworkers. In Sect. Ill, we define the disordered DTFB model, and analyze some realizations 
of disorder at fixed strain. The experimentally more relevant case of fixed stress is discussed 
in Sect. IV. Section V is devoted to the discussion of the Edwards- Anderson parameter, 
whereas some final comments and a summary of the results are given in Section VI. 



II. THE TFB MODEL 



The homogeneous Thermodynamic Fiber Bundle (TFB) model consists of N fibers (har- 
monic springs) subjected to the same load. Let the set of variables {tf^i — 1, ...N} describe 
a configuration with broken (ti = 0) and intact {U — 1) fibers. The energy (Hamiltonian) of 
a configuration of this system is given by 

^ N N 

H^-ke'Y^U-DY^U, (1) 

i=l 1=1 

where /c > is an elastic modulus, e is the common strain, and D > is the (uniform) 
dissociation energy. We assume that the breaking process is reversible, with probabilities 
related to the Boltzmann factors, and write a canonical partition function in the "strain" 
ensemble. 



Z = ^exp 
{ti} 



^ N N 



i=l 



1 + exp ( --(3ke^ + (5D 



i=l 



(2) 



where /9 = 1/T is the inverse of the temperature. From this partition function, we have the 
Helmholtz free energy per fiber. 



1 + exp ( ~(3ke'^ + (3D 



(3) 
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Given a configuration {tj}, the fraction of intact fibers is written as 



1 



(4) 



whicli leads to tlie canonical average 

(0)t = 



i=l 



1 



(5) 



1 + exp {f3key2 - [3D) 

as in the original work of Selinger et al. 

In order to make contact with the calculations of Selinger and coworkers, we can rewrite 
the canonical partition function as an integral over the fraction of intact fibers 0, 

'-^ke'<P + D^ I (6) 



{u} \ i=i J ^ 



Using an integral representation for the Dirac 5-function, and discarding irrelevant terms in 
the thermodynamic limit (large iV), it is easy to see that 



where 



.exp [-(3NT] 



r={he^-D)cP+^[{l-<P) ln(l-0) + 01n, 



(7) 



(8) 



The minimization of JF with respect to (j) leads to the Helmholtz free energy of Eq. ([3]). 

We now consider the more realistic situation of fixed stress a. In the stress ensemble, we 
write the partition function 

N 



Y = de exp {-(3aeN) Z = ^exp 

{ti} 



i=l 



Again, if we use an integral representation of the (5-function, and discard irrelevant terms, 
it is easy to write 

Y= f d(i)exp[-pNg], 



(9) 



with 



2 -1 

-^-D0+-[(l-0)ln(l-0) + 01n. 



(10) 



whose minimum with respect to leads to the expression of the Gibbs free energy per fiber 
bond. Sehnger and coworkers point out that the function Q in terms of 0, for sufficiently 
small values of the stress cr, displays a local relative minimum at 0*, with < 0* < 1, which 
corresponds to a mctastable state with a fraction 0* of intact fibers. Also, it is shown that 
increasing the stress a, there is a decrease of both the values of 0* and of the height of the 
barrier separating the local minimum from the global minimum of the free energy at = 0, 
with a fracture threshold at the spinodal (Tc = (7c (T, -D, k). 



III. THE DISORDERED DTFB MODEL 



We now assume that each fiber bond is characterized by a different dissociation energy 
Di, for i = 1, ...N. The disorder is fixed, quenched, so that {Di} is a given set of independent 
and identically distributed, random variables, with a probability distribution p{Di). The 
energy of a configuration of this system is given by 



N 



N 



H=^ke^J2^i-J2^iti- (11) 

1=1 i=l 

Again, we assume a thermodynamic process, at temperature T, with Boltzmann weights. 
In the canonical strain ensemble, given the configuration {Di}, we write the partition func- 
tion 



N 



^{A}=n 



1 + exp ( --pke^ + pD, 



(12) 



In the thermodynamic limit, we can use the law of large numbers in order to write the 
(self-averaging) Helmholtz free energy 



1 + exp ( --pke'' + pD 



where 

(•••) = J dDp{D) 
The average fraction of integer fibers is given by 

0(6, T) ' ^ 



(13) 
(14) 

(15) 



1 + exp {l/3ke^ -PD) /' 
which resembles the form of the density of states of free fermions. In the zero temperature 
limit, this fraction becomes a step function, so that 

1 



lim0(e,T) = l-P(-A;6^), 



(16) 



where P{D) is the cumulative probabihty distribution, 



D 



P{D) ^ / dD'piD'), (17) 
Jo 

in agreement with the predictions for the irreversible FB model at constant strain. 

In the following subsections, we consider a situation of fixed strain, with uniform and 
bimodal distributions of the dissociation parameter D. 



A. Uniform disorder 

In the usual version of the FB model {2I, one assumes a uniform distribution of dissociation 
energies and identical harmonic elastic terms. We then compare the main features of the 
disordered DTFB model with the corresponding results for the uniform TFB and the FB 
models. In the strain ensemble, it is straightforward to calculate the averages given by Eqs. 
f|T3|) to f|T7|) . In particular, we look at the average fraction of intact fibers, and the stress- 
strain (cr versus e) characteristic curve. If the dissociation energies are uniformly distributed 
over the interval [a, b], with < a < 6, we have 

^1 + g-(|fce2-b)/r" 



^ 1 

— a 



1 + e 



-ake'^-a)/T 



(18) 



In Figs, la) and lb), we plot as a function of temperature for different values of the strain. 
The corresponding behavior of the uniform TFB model is shown in Figs. Ic) and Id). In 
Fig. (2), we plot as a function of strain, e, for different values of temperature T. The 
main differences are at low temperatures. In this regime, the presence of disorder smooths 
the transition from = 1, at low strain, to = 0, at large strain. At T = the transition is 
abrupt in the absence of disorder, but it becomes continuous in the disordered DTFB model. 
In this case the curve is directly related to the probability distribution for the dissociation 
energies, Eq. (1161) . and reproduces the behavior of the usual FB model, with a discontinuity 
in the derivative. At low but non vanishing temperatures, there are no singularities in both 
models, but in the TFB model drops more quickly with increasing strain. This is reflected 
in the stress-strain characteristic curve a = kipe shown in Fig. 3a). 
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B. Bimodal distribution of disorder 



Although the case of uniform disorder does not differ too much from the case with no 
disorder, there are strong differences with other forms of the distribution of disorder. For 
example, we may consider a bimodal distribution, 



p{D) = (1 - c)6{D - Di) + c6{D - D^). 



(19) 



where c, Di, and D2 > Di are non negative parameters. In the strain ensemble, the 
Helmholtz free energy is thus given by 



/ 



1 

1 



'1 - c) In 



1 + exp — /3ke^ + (3D] 



icln 



1 + exp ( --I3ke^ + I3D2 



with the average fraction of intact fibers 

1 - c 



(20) 



(21) 



^ 1 + exp (i/?te2 - pDi) 1 + exp [\l3ke'^ - /3D2) ' 

In Figs. 4a) and 4b), we draw graphs of versus temperature T for different values of the 
strain, with c = 0.7, Di = and D2 = 1. In contrast to the case of uniform disorder, 
the graph of the average fraction (p displays pronounced extremal points, which give rise to 
a counter-intuitive behavior. At small values of the strain, the average fraction of integer 
fibers first increases with temperature, and then drops slowly. At large values of the strain. 
Fig. 4b), it first decreases and then increases monotonically at higher temperatures. This 
behavior can be understood if we consider that in this case behaves as a linear combination 
of two TFBs with no disorder and different dissociation energies. In each TFB is increasing 
for the strain e less than a certain threshold Eq, and decreasing for e > Eq, with Eq = 

1 /2 

{2D/k) ' , Eq. With bimodal distribution, we anticipate that should display a more 

1 /2 

or less pronounced extremum for Ei < E < E2 (with Ei = {2Di/k) ), since in this region 
there is a superposition of monotone increasing and decreasing curves. In fact, the form of 
the disorder distribution is crucial for determining the mechanical properties of the system. 
In particular, characteristic load curves may differ widely, depending on whether the disorder 
distribution has a finite support, short tails, long tails, and other geometric features. A study 
of this dependence will be the subject of a forthcoming paper [22|. 
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IV. THE DTFB MODEL IN THE STRESS ENSEMBLE 



In most of the real experimental situations, the external stress a is given, instead of 
the internal strain e. As in Section III, given a disorder configuration {-Di}, we write the 
partition function in the stress ensemble. 



Y{D,} = j dee^-^{-paeN)Z{Di} 

N 



^exp 



i=l 2/cX]j=i^i 
In order to carry out the calculations it is useful to introduce a new variable 

N 



i=l 



(22) 



(23) 



which is a typical variable for disordered systems, representing a sort of disorder-averaged 
magnetization, the projection of the system configuration ti upon the disorder. Using Q and 
the definition of 0, given by equation (jl]), we write 



exp <^ (3N 



Q 



a 



2k(t) 



(24) 



We now introduce integral representations of the 5-functions. In the thermodynamic limit, 
we invoke the law of large numbers in order to see that the final expression for the free 
energy is self-averaging with respect to disorder. Thus, we have 



with 



where 



Y {Di} = Y = j d(p J dQ exp [-PNG] 

g = -Q 



a' 1 

2k(f) ~ (3N 



Inl], 



Vt = j j dxdyexp [NJ^(x,y)], 
T (x, y) = 0x + + (In [1 + exp (— x — yD)\) 



(25) 
(26) 

(27) 
(28) 



and (...) indicates an average with respect to the probability distribution p{D). A saddle- 
point calculation leads to an asymptotic expression of Q, which in turn leads to the final 
(Gibbs) free energy of this system. For simple forms of the probability distribution we can 
even write some analytic expressions. 
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A. Bimodal distribution of disorder 



For a bimodal distribution, given by Eq. (fT9|) . we can write analytical expressions. In 
particular, with Di = and D2 = D > 0, and defining 

X = §, (29) 

we have 

g = g{T, a; x) = -Dx-^- (0, x) , (30) 

where 

s (0, x) = - [(1 - c) + X - 0] In [(1 - c) + X - 0] + (X - 0) ln(x - 0) + 

+ Xlnx - (c- x)ln(c- x) + cine + (1 - c) ln(l - c). (31) 

In Fig. [5l we draw equilibrium values of x and for c = 0.7 and applied stress a = 0.3. It is 
interesting to see that there are two distinct values of x for each value of 0, corresponding 
to different temperatures. This implies that there are two different microscopic states at 
each point of the load curve (cr, e). The turning point around T = 0.1 corresponds to 
the extremum of in Fig. [61 Figure [7] reports the dependence of on o" for different 
temperatures. In the next Section we make some additional comments on these results. 



B. Uniform disorder 



Let us look at the case of a uniform distribution of dissociation energies in the interval 
[a,b]. Since we can no longer express the free energy in terms of elementary functions, 
we are forced to resort to some numerical calculations. In Fig. [8], we compare the phase 
diagrams of the disordered DTFB and the uniform TFB models. The solid line represents 
the spinodal for the Helmholtz free energy. In the DTFB model, the averages are taken with 
respect to two different uniform distributions. It is seen that the main effect of disorder 
is the weakening of the system at low temperature. In the case of wider distributions of 
disorder, the spinodal displays a plateau as T — > 0, implying a constant failure stress in that 
region. For temperature larger than 1, the TFB and the DTFB models display identical 
behavior, and the disorder does not affect the system. 
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V. THE EDWARDS-ANDERSON PARAMETER 



In the absence of disorder we may take as the only order parameter of the system. We 
may replace the ensemble average of a quantity {A) by a time average lim^-^oo / A(t)dt, 
and thus interpret the average fraction of intact fibers as the average fraction of time during 
which a fiber remains intact, corresponding to the thermal average (ti)T- This is the same 
for all the fibers since all of them are identical. If disorder is introduced, and the dissociation 
energy of each fiber becomes different, this no longer holds. The thermal average may be 
different for each fiber, depending on the value of D. Thus, the same macroscopic state, (j), 
may correspond to different microscopic states, and it might be interesting to distinguish 
among them, even for practical purposes. A quantity which can operate this distinction is 
an analog of the Edwards- Anderson parameter |23| , 



At constant strain, we can write 



q{e,T) = dDp{D) ^, (32) 



from which it follows that 



lim q = lim di^ = -. 



2 



as it could have been anticipated. Then, it seems useful to introduce the quantity r = q 
as a measure of the degree of "freezing" of the system. 

Figure [9] shows a detail of the strain stress curve of Fig. [3] for the bimodal distribution 
at different temperatures. This plot shows that the curves for T = 0.1 and T = 0.3 cross at 
a point in the {e, a) plane. This same macroscopic state corresponds to different underlying 
microscopic states characterized by different values of the elastic modulus c. These values 
can be related to the microscopic state by the analog of the Edwards-Anderson parameter. 
In fact we have: 

c=^ = k4>~k'-{<p-q). (33) 

From this expression, together with a = k(j)e, we obtain the values of and q. In particular, 
at the failure, c = 0, we have 

(34) 



^ ke 



T 

ke^ 



10 



Another instance of this lack of uniqueness is given in Figs. [TO] a) and b), which have 
been generated by simulated annealing. The graphs in these figures show the dependence 
of the "degree of freezing" r = g — 0^ on the applied stress at some low temperature values. 
The corresponding curves for cj) are given in Figs. [7] a) and b). In analogy with the stress vs 
strain diagram, the curves for different temperatures cross at some points, where however r 
assumes distinct values. 



VI. CONCLUSIONS 



We have shown that the introduction of disorder in the Thermodynamical Fiber Bundle 
model 2^ affects its behavior in many respects. Besides the onset of some expected features, 
as a decrease in the failure stress at low temperatures, there appear some new features 
associated with the distribution of disorder. For a bimodal disorder distribution, which 
leads to an analytically tractable problem, the fraction of integer fibers is non-monotonic in 
terms of temperature, with an extremum at low temperatures. As a consequence, at either 
constant stress or constant strain, the system may display the same values of the fraction of 
integer fibers at different temperatures. However, the corresponding microscopic states are 
different, which is relevant for understanding delayed fractures, and can be characterized by 
different values of an analog of the Edwards Anderson parameter. 

A.V. is grateful to Fergal Dalton for many useful discussions and suggestions. A. P. S.R.S 
acknowledge financial support from CNPq and CNR. 
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FIG. 1: Fraction of integer fibers ^ as a function of temperature for different values of the applied 
strain e: a), b) in the DTFB model with a uniform distribution of disorder in [0,2]; c), d) in the 
homogeneous TFB model (with no disorder) with D = 1. For 0.9 k, e k, \f2 these curves differ 
very little from each other. We always take k = \. 
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(a) (b) 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



8 e 

FIG. 2: Fraction of integer fibers as a function of strain at different temperatures: a) in the DTFB 
model with a homogeneous distribution of disorder with hmits and 2; b) in the homogeneous 
TFB model (no disorder) with D = 1. We always take k = 1. 




FIG. 3: Stress-strain characteristic curves at different temperatures: a) disordered DTFB model 
with a homogeneous distribution of disorder in [0,2]; b) homogeneous TFB model (no disorder) 
with D = 1; the case T = corresponds to the ordinary FB model; c) DTFB model for bimodal 
disorder, Eq. p^ . with parameters c = 0.7, Di = and D2 = 1. We always take k = 1. 
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FIG. 4: Fraction of integer fibers (p versus temperature T for a bimodal distribution, Eq. ()19p . with 
Di = 0, D2 = 1, c = 0.7 at constant strain with k = 1. 
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FIG. 5: Parameters x vs (p, after minimization of Eq. (j30p and from numerical simulations, in the 
case of bimodal disorder, Eq. (jl9p . with c = 0.7 and stress a = 0.3, from Eq. (j30p . and k = 1. It is 
seen that x assumes two different values for the same value of (p. 
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FIG. 7: Stress dependence of cp at various temperatures for the bimodal distribution of disorder 
with Di = 0, D2 = I, c = 0.7, and a Young modulus k = 1. 
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FIG. 8: Comparison between the phase diagram of a case without disorder (TFB with D=l) and 

the phase diagrams for two homogeneously disordered cases: Di G [0.5; 1.5] and G [0;2]. The 
Young modulus is k = 1. 
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FIG. 9: A bimodal distribution of disorder with Di = 0, D2 = 1, c = 0.7 (we always take k = 1). 
It is seen that for two different temperatures, T = 0.1 and T = 0.3, stress and strain assume the 
same values (around e = 0.35 and a = 0.28). The different microscopic configurations can be 
distinguished by different values of the Edwards- Anderson parameter q. 
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FIG. 10: Stress dependence of r at various temperatures for the bimodal distribution of disorder 
with Di = 0, D2 = 1, c = 0.7, and a Young modulus A; = 1. The inset in a) shows that curves for 
distinct temperatures have distinct values and do not intersect. 
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